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between the Hamiltonian approach developed here and the known variational one is 
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1. Introduction 

Riemannian cubic polynomials (RCP), also called Riemannian cubics, can be seen as a 
generalization of cubic polynomials in Euclidean spaces to Riemannian manifolds. The 
cubic polynomials on a Riemannian manifold are the smooth solutions of the fourth 
order differential equation 

+ Rhn^'37 37 = 0, (1) 



dt^ V dt^ ' dt y dt 

where D/dt denotes the covariant differentiation and R the curvature tensor. The 
equation ([1]) is the Euler-Lagrange equation of a second order variational problem with 
Lagrangian given by | (D^x/dt^, D^x/dt^), where (-, ■) denotes the Riemannian metric. 
This variational problem was first introduced in 1989 (see [23]) and explored from a 
dynamical interpolation perspective in 1995 (see [15]). Interesting points related to this 
subject have been developed in the last few years, namely a geometric theory surprisingly 
close to the Riemannian theory of geodesies (see[2l[3lll[T0l[l2l[T3l[Tl[T7l[26l|271[2g 
We recall, in particular, a result which says that if V denotes the velocity vector field 
of a cubic polynomial x, then 

1 /BV BV\ /DV ^ \ 



2 \ dt ' dt / \ dt2 
is invariant along x. In Riemannian context, Ji plays a role similar to the one played by 
the length of the velocity vector field in the theory of geodesies (see, for example |14j). 
Recently, in [3], [26], ET] |29] , the analysis of RCP from a variational point of view was 
carried out for locally symmetric manifolds and a second invariant was obtained: 

B^V DV\ /DV BV\ 

The analysis of RCP given in [3l [26] is qualitative, with special attention to the case of 
the Lie group 5*0(3), where RCP correspond to Lie quadratics on the Lie algebra. The 
article [3] introduces a reduction of the RCP's equation for this Lie group of rotations. In 
[26] some results on asymptotics and symmetries of cubics are proved for the particular 
case of the so-called null cubics on 5*0(3). In [27], the author solves by quadratures the 
linking equation on 50(3) and 50(1,2) of the Riemannian cubics. Finally, [29] studies 
n-th order generalizations of RCP introduced in [12] . 

To our knowledge, the first Hamiltonian description of the RCP problem has been 
considered in [13] (made in collaboration with one of the authors). The present paper 
deals, for the case of arbitrary compact and connected Lie groups, with a different 
Hamiltonian description of the problem. Here we use a presymplectic approach to the 
Pontryagin's Maximum Principle inspired by some ideas of [TJ [9l [161 122]. Namely, 
we consider the intrinsic geometric approach used in [T6| [22] for a first order general 
optimal control problem, and similarly considered in [7] for time-dependent optimal 
control problems by using the jet bundles framework. In a similar way, [9] gives the 
geometric treatment of the Lagrangian dynamics with higher-order constraints. The 
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description of RCP (on an arbitrary manifold) using these geometric ideas were first 
presented in |1] by the authors of the present paper. Recently, in [5], the authors have 
treated the particular situation of the dynamic control of the spherical free rigid body, 
a mechanical system with configuration manifold given by the Lie group SO (3). The 
new contribution of this presymplectic formalism is to use the Lie group structure of the 
semidirect product of the Lie group and its Lie algebra, Gxg, which is the state space of 
the optimal control problem. This allows us to use classical results from [H |8], adapted 
to this Lie group structure. Namely, we present the left trivialization of the Hamiltonian 
system, a set of equations which lives in the manifold given by the cartesian product of 
the semidirect product above mentioned and the dual of its Lie algebra, G x gx g* x g*. 

The main goal of this work is to reduce the degrees of freedom of the left trivialized 
Hamiltonian system. We first apply the symplectic point reduction theorem (|23l |28] ) 
and then explore the reduced dynamics using a suitable symplectomorphism. The 
reduced Hamiltonian vector field lives in the manifold given by the cartesian product of a 
coadjoint orbit, the Lie algebra and its dual, Orj x g x g* . Furthermore, some invariants 
along the extremal trajectories are characterized as a crucial point to develop, in a 
future work, a study of the integrability of the Hamiltonian system. In fact, using the 
Lie-Cartan theorem ([6]), we obtain an interesting result on the number of independent 
integrals of motion in involution. 

The plan of the paper is as follows. Section 2 is a support of the other sections 
and it reports some notes on compact Lie groups and fixes the notation. Section 3 
begins with the introduction of the optimal control problem of cubic polynomials and 
presents the corresponding presymplectic approach. After that, we provide the left 
trivialized description equivalent to the variational one ([I5]) in a similar way to what 
happens in [13]. However, it is important to remark that our Hamiltonian system and 
the one in [13] are different. The last section contains the reduction of the Hamiltonian 
problem by means of the Marsden-Weinstein technique. As we have mentioned above, 
the system will be reduced to a system on the cartesian product of a coadjoint orbit, 
the algebra of the Lie group and the dual of this algebra. In this context, besides the 
geometric deduction of the two known invariants, we find some more invariants along 
the Riemannian cubics and find a relevant result on its independence and involution. 

2. Preliminary notes 

The present section gives some notations, definitions and results from the Lie groups 
theory, which we will be using throughout the paper. 

2.1. The Lie group G 

Let G be a connected and compact Lie group with identity e. The corresponding Lie 
algebra is (g, [.,.]), where [.,.] is the bracket operation and g* denotes the dual space 
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of this algebra. Furthermore, the elements of G are denoted by x or (7 and the maps 
G X G ^ G,{x,g) (— )■ xg and G ^ G,x x~^ are the multiplication and inversion 
operations for the Lie group G, respectively. Given x,g E G, let : G G and 
Rx '■ G ^ G be, respectively, the left and right translations by x. The tangent of at 
g is denoted by TgLx and T*Lx represents its transpose. Recall the following definitions: 

- The adjoint representation of the Lie group G is denoted by Ad. It gives for each 
X G G* an algebra automorphism defined by Ad^ = Te(-R^.-i o L^). 

- The adjoint representation of the Lie algebra g is the tangent of Ad at the identity 
e and it is denoted by ad. For each F, Z G g, we have adyZ = [F, Z]. 

- The map Ad* : G ^ Aut{Q*) defined, by [Ad*(a;)](0 := Ad^-i^ = { o Ad^.-i, for 
each X E G and G 0*, is called the co-adjoint representation of G. 

- The co-adjoint representation of Q is the map ad* : g — )■ Aut{Q*) defined, for each 
y G g and e e g*, by [ad*(r)](0 := -ad^^ = -i o ady. 

Since the Lie group is assumed to be connected and compact, we can guarantee the 
existence of a bi-invariant metric on G, which we shall denote by (.,.). This statement 
and the following result can be found for instance in |24j . 

Theorem 1 [24^ If G is a Lie group equipped with a bi-invariant metric, the metric 
connection V and the curvature tensor R associated with that metric are given, 
respectively, by VyZ = ^[Y, Z] and R{Y, Z)W = -l[[Y,Z],W], where Y, Z and 
W are left invariant vector fields. Furthermore, the first above equality implies that 
{[Y,Z],W) = {Y,[Z,W]). 



In the course of this paper, we shall fix an orthonormal basis in the Lie algebra g. 
The corresponding dual basis is a basis of the dual space g*. These two bases generate 
left invariant frame and co-frame fields on G, respectively. We assume the following 
notations: 

- Let y be a curve in g and ^ a curve in g*. We represent by Y (respectively, 
^) the element of g (respectively, g*) which has components with respect to the 
left-invariant frame (respectively, co-frame) field above mentioned, given by the 
derivative of the components of Y (respectively, ^). 

- Given ^ G g*, the tangent vector identified with this co- vector by the Riemannian 
metric will be denoted by G g. That is, ^(Y) = (X^, y), VF G g. 

- With the above notation, it is simple to verify that = X^ and X^.^^^ = —adyX^. 
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2.2. The tangent bundle Lie group and its left trivialization 

Lemma 1 ^ilH]/ The tangent bundle TG is a Lie group with a group operation defined as 
the tangent prolongation of the original one on G. That is, the multiplication operation 
for TG is defined by 

(Vx, Vg) G T^G X TgG I )■ VxVg TxRgVx ~\~ TgLxVg G TxgG 

and the inversion is defined by 

Vx G TxG I — yv~^ = - {T^Lx-i o TxRx-i) (vx) G T^-iG. 

Consider the semidirect product G(s)g of the Lie group G and the Lie algebra g 
regarded as abehan group, under the right action of G on g, (x, F) G G x g i — )■ Adx-iY. 

Lemma 2 The semidirect product G(S)g is a Lie group whose underlying manifold is 
the cartesian product G x q and group multiplication law 

{x,Y){g,Z) = {xg,Adg-iY + Z), 

for [x, Y), {g, Z) G G x g. The inversion is defined as {x, Y)~^ = {x"^, —AdxY). 

In order to simphfy the notation, we shall denote the semidirect product by G x g 
instead of G(s)g. 

The semidirect product structure considered here is a special case of the general one 
defined by a right representation of a Lie group on a vector space that may be found in 
works on semidirect products, particulary the ones on models of continuum mechanics 
and plasmas where it is convenient to work with right instead of left representations 
(see, for example [20]). 

Proposition 1 f21^ The left trivialization ofTG determined by the map 
X: TG — y Gxg 

Vx^TxG I — > {x,TxLx^iVx) , 
allows us to write the Lie group diffeomorphism TG ^ G x g. 

We introduce now some important notations used in the rest of the paper. 

- The elements of the tangent bundle T(G x g) are denoted by 

{vx, y, U) G T(,,y)(G X g) = TxG x {Y] x g. 

- The second tangent bundle of G, T^G, can also be trivialized by using the map A 
and then realized as a vector bundle over Gxg, which is a vector subbundle of 
T(G X g). We represent this bundle by T^G and denote its elements as 

iyx-, U) G T(^x,TxL^-iVx) 

(G X g) ~ TxG X g. 

- The elements of cotangent bundle T*{G x g) are represented by 

G T(;,y)(G X g) = t;g x {Y} x g*. 



(4) 
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In the previous statements we are considering (x, Y) e G x g. Throughout this paper we 
will, for the sake of simplicity, occasionally assume the identification between elements 
of T(a;y)(G X g) (respectively, T*^y){C^ ^ s)) elements of T^G x g (respectively, 
t;g X g*). 

According the Lie group structure chosen in lemma [2], we easily compute: 

T(e,o)L(x,Y)iZ, U) = (TeL^Z, U + adyZ) (5) 

and 

ad(y,z) (Y', Z') = (ady r, ady Z' + adzY') , (6) 

where x G G and Y, Z, U, Y', Z' G g. Obviously, these formulas can be derived from the 
general known ones from the theory of semidirect products. 



3. Hamiltonian system 

The aim of this section is to give a Hamiltonian description of the optimal control 
problem of cubic polynomials on G based on some material published in [1], where 
we used a geometric formulation similar to the one developed in [9] for higher order 
constrained variational problems. The section begins with the introduction of the 
optimal control problem, where by means of the left translation on G, the state 
space has been left trivialized to be G x g instead of TG. After that, we apply a 
presymplectic constraint algorithm and the result is a Hamiltonian system on a space 
symplectomorphic to T*(G x g). Using again a left trivialization, but now determined by 
left translation on the group Gxg, we pass to a Hamiltonian description on Gxgxg* xg*. 



3.1. Optimal control problem 

Considering the left trivialization (jl]) of TG, the state space for our problem may be 
taken to be the semidirect product Gxg and the bundle of controls as the second 
tangent bundle T^G. The optimal control problem of cubic polynomials on G consists 
in finding the G^ piecewise smooth curve 7 : [0, T] — )■ T^G with fixed endpoints in state 
space, mimmizmg the functional /p L(7(t))dt, with T G M+ fixed, for L : T^G M the 
cost functional defined by 

L{v.,U) = \{U,U) (7) 
and satisfying the control system 

A(r,^(^(t))=F(7(t)), (8) 

where : T^G — t- G x g is the natural projection and F : T'^G — > T(G x g) is the vector 
field along this projection defined by 

U). (9) 
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T^G ^T(Gx0) 




Notice that, according the notation set in the subsection 12.21 a curve 7 in T'^G is 
defined by means of three elements: a curve x in G; a vector field along x (which 
can be seen as a curve in TG satisfying vtc o K^. = x, where ttg '■ TG — )■ G is the 
canonical projection); and a curve U in g. So, 7(t) = (^^(t), f/(t)) G Tx(t)G x q and 
we have rli^it)) = {x{t) ,Tx(t)Lx{t)Y:c(t)) . Consequently, using the appropriated basis of 
left invariant vector field on G to develop the calculus, it is simply to prove that the 
control system ([8]) can be written as 

BY 

x{t) = n(t) -^{t) = nL,^t)U{t), (10) 
which is a version of the control system presented in papers fLli [T3] . 



3.2. Dynamics of the optimal control problem 

The co-state space of our system is the cotangent bundle T*{G x q). The dynamics of 
the control problem is described by a presymplectic system (T, fl,lTj whose total space 
is the vector bundle over G x q given by 

r = T*{GxQ)x^^^T^G. (11) 

The elements of this space are points in T*G x {Y} x q* x g denoted by {ax,Y,C,,U), 
where {x,Y) E G x q. Consider the canonical projections pri : T ^ T*{G x g), 
(a,,F,e,f/) ^ («x,l^,0 and pr^ : T ^ T^G, {a,,Y,i,U) ^ {T,L,Y,U). The closed 
two form is defined by the pull-back 

Q = {pr,yQ,, (12) 

with Qi denoting the canonical symplectic two form on the space T*{G x q). The 
Hamiltonian is defined by H =-C pri,F o pr2 ^ —L o pr2, where F and L are defined 
by (E]) and Qj and ., . ^ stands for the canonical duality product of vectors and 
covectors on G x g. Then, 

H{a,, Y, e, U) = {T:L,a,) {Y) +i{U)- ^(t/, U). (13) 

The dynamical vector field of the system is the vector field Xjj : T — t- TT solution of 
the dynamical system ixj^^ = dH. 

Notice that the optimal control problem is obviously regular and thus applying the 
geometric algorithm of presymplectic systems (see [HI [19]) to (T,^l,H), we obtain a 
symplectic system on the manifold Wi = {{ax,Y,^,U) E T : U = X^}, where E g 
is the tangent vector identified with the covector E g* hj the Riemannian metric 
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of G. Hence, (VTi, i^vvi) is a symplectic system, with and H^-^ being the 
restrictions to Wi of f|T2l) and f|T3l) . respectively. The map / defined below, gives us a 
diffeomorphism between the symplectic manifolds (T*(G x and 

14 

So, we have a symplectomorphism between the two manifolds (see [1], p 177). In this 
sense, we construct the Hamiltonian Hi := Hwi o / : T* (G x g) — )■ M. We get 

Hi{a,,Y,0 = {T:L,a,) (F) + ^^^i) , (15) 

for each ^) G T^^y){C^ ^ s)) where {x,Y) E G x q. Furthermore, the existence 

of the symplectomorphism f|T^ allows us to conclude that (see [1], p 194) the study of 
the dynamical system defining the Hamiltonian vector field Xjj^ associated to H^-^^ is 
reduced to the study of the system 

tx„^ni = dHi, (16) 

where the vector field Xh^ : T*{G x q) ^ T{T*{G x g)) is the push- forward of Xjj^. by 
f~^, = (/~^)*X-^^^. . The integral curves of this vector field determine, when lifted 
back to the original manifold, the solutions of the optimal control problem (see [T6l [7]). 



3. 3. Left trivialization of the dynamics 

Consider now the left trivialization of the cotangent bundle T* {G x q), determined by 
the diffeomorphism p defined from T*(G x g) to the space G x g x g* x g* as 

P(«x,^,0 = (a;,>^,r(e,o)^(x,y)(«x,^,0) , 
which using ([5]) gives p{ax,Y,^) = {x,Y,T*Lxax + adyC,,^), for each (a^,y, ^) in 
^(x,y)(^ X g), where (x, Y) e G x q. Observe that, if (x, F, /i, ^) G G x g x g* x g*, then 
p"i(x, y, /i, = (T;L,-i (/i - ad;.o , e Tl^^Y)iG x g). 

The left trivialization of the Hamiltonian (fT5|) is given by H := Hio p~^. We easily 
conclude that 

H{x,Y,p,0 = KY) + laX^). (17) 

Since p is a diffeomorphism, we can endow (see pQ, p 177) G x g x g* x g* with a 
symplectic structure, as = {p~^)*Qi. Furthermore, p is a symplectomorphism and 
(see [1], p 194) the Hamiltonian vector field X^i defined by (fT6l) may be left trivialized 
to G X g X g* X g* by considering the push-forward by p of the Hamiltonian vector field 
associated to Hi, Xh ■= p^Xu^. 

The proposition below leads us to the expression of Xh: 

Xh{x, r, p, = (TeL.r, x^, 0, -p + ad;.o, (is) 

for each (x, Y, p,S,) E G x q x q* x g*. 
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Proposition 2 The following set of differential equations describe the motions of the 
Hamiltonian system (G x g x g* x g*, f2, i7) 

/i = U 

^ = — /i + ady^. 

Proof. Let z = {x,Y,fi,^) be an integral curve of Xh- Following ([8], section A. 3, 
example 3) the Hamiltonian's equations on the cotangent bundle of the Lie group G x q 
are 

•\ ^ , fdH ^ , dH ^ 

X,Yj = T(e,o)L(,,y) ( ^W'^l^ 

In this notation, dH{z)/dx is regarded as an element of T*G, dH{z)/dY as an element 
of 5*, dH{z)/dfi and dH{z)/d^ as elements of 0. Use ([5]) and (12 .2^ to rewrite the 
previous system as 

i = TeLx^{z) 

Y = f{z)+ady^{z) 

fi = -T:L^mz) - ad;.|f (;a;) + a.d*^ + a.d*aH,.^ 

From the expression of the Hamiltonian function ( fT7|) it comes dH{z)/dx = 0, 
dH{z)/dY = /i, dH{z)/dji = Y and dH{z)/d^ = X^. Now, substitute these expressions 
in the above system, use the fact that ad^^^ = and the result follows. □ 

Remark 1 It will now be interesting to see how the dynamics described by l[T8\) is related 
with the known variational approach of cubic polynomials. To proceed, we begin with the 
following remarks: 

• First write the last equation of [T^) as an equation on the Lie algebra, using the 
identification of covectors and tangent vectors giving by the Riemannian metric ofG 
(see the end of subsection \2. 1\ for details on notation). We get X^ = —X^ — adyX^. 

• Differentiate the above equation and use the third equation of / fTPj) . to obtain 
X^ + adyX^ + adyX^ = 0. Use the second equation of [W) . to get y F] = 0. 

We have just shown that each solution of the equations of Hamilton ( figj) gives rise to a 
solution of the equations 

Y = T.,L,-ix Y+[Y,Y] = Q. (20) 

Conversely, solutions of ^2Wl satisfying Y = X^ and X^ + X^ + adyX^ = 0, correspond 
to solutions of / flgj) . 
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The equations are the Euler- Lagrange equations (J\) that define the cubic 

polynomials on a Lie group, which were proved in ITB^ as an extension of the proof 
that had already been given in f25^ for 5*0(3). (The proof use some facts derived from 
theorem [II ) 

4. Reduction of the Hamiltonian system 

The purpose of this section is to study the symmetries of the Hamiltonian system 
{G X Q X g* X g*,Q,H) described in the previous section and use that to reduce the 
corresponding dynamics, ehminating degrees of freedom in the system. The idea is to 
apply the symplectic point reduction theorem (see |23] for the original references and 
[28] for full details in this subject) and carry out the appropriate interpretation of the 
reduced Hamiltonian system for the study of important questions as the integrability 
of the system. Namely, we shall focus our attention on the integrals of motion of the 
reduced Hamiltonian system. 

4-1- Symplectic point reduced space 

Let (p be the smooth left action of the Lie group G on G x g x g* x g* defined by 

0((7,(x,F,/i,O) = (^7a;,r,/i,O, (21) 

for each g & G and (x,y, /i,^) & G x g x g* x g* . The moment map of is the map 
J : G X g X g* X g* ^ g* defined, for each (x, F, ^) G G x g x g* x g*, by 

J(x,y,/i,0 = Ad^i(^-ad^O- (22) 

The action cf) can be seen as the left trivialization to G x g x g* x g* of the cotangent 
lift of the action of G on G x g given, for each g E G and (x, Y) G G x g, by 
{g, (x, Y)) L(^g^Q){x, Y) = {gx, Y). Recall that every cotangent lift action is symplectic 
and has momentum map Ad*-equivariant (see [1], p 283). So, it is easy to verify the 
following statement: 

(A) is a symplectic action with momentum map Ad*-equivariant. 

Observe now that the action </> is proper since it is an action of a compact Lie group. 
Moreover, is obviously free and hence the symmetry algebra of every point in 
G X g X g* X g* is zero, which is equivalent to saying that every ?7 G g* is a regular 
value of the momentum map J. 

Let ?7 G g*. Consider the coadjoint isotropy subgroup of rj, defined by 

G,:={geG : Ad*-ir/ = r/} (23) 

and also the level set J~^(?7) of the momentum map J. Note that 

J-\r]) = {(x, F, /i, G G X g X g* X g* : /i = Adlr] + ad^O- (24) 

Because is a symplectic G-action on the symplectic manifold G x g x g* x g* and 
?7 G g* is a regular value of J, we see that J~^{rj) is a submanifold of G x g x g* x g*. 
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Furthermore, as a consequence of J being Ad*-equivariant, we easily prove that J~^{ri) 
is G^-invariant. The comments now exposed allow us to conclude that acts on J~^{ri) 
and that the orbit space 

(G X X r X 0*)^ := J-\v)/G, (25) 

is well defined. The action of G^ on J^^{ri) is obtain by restriction of to the subgroup 
(l23l) and to the G^-invariant submanifold (Elj). It turns out that the action </> is proper 
and free and that by definition G^ is a closed subgroup of G, thus (see [28], p 60): 

(B) The action of G,, on J~^{ri) is proper and free. 

This result guarantees that the orbit space ( 125|) is a smooth manifold and that the 
corresponding projection map is a surjective submersion. 

Since the conditions (A) and (B) are satisfied, we are able to apply the symplectic 
point reduction theorem. The theorem states the following result: 

The reduced space (G x g x g* x g*)^ has a unique symplectic structure fi^ 
characterized by the identity vr*f2^ = z*^, where is the canonical inclusion 
from J^^{ri) to G x q x g* x g* and tt^ is the projection of J^^irf) onto 
{Gxgxg* X g*)^ . 

The symplectic manifold ^(G x g x g* x g*)^, Vl^^ is called the symplectic point reduced 
space at 77. 



Let us now explore in more detail the reduction obtained. More specifically, we 
will interpret the symplectic point reduced space in a strategic way to conduct further 
studies. In what follows, we shall adopt the notation for the above G^-action on 
J~^{r]). First, we notice that from (IMl) . the submanifold J~^{ri) is diffeomorphic to the 
semidirect product Gxgxg* (of the Lie group G x g and the vectorial space g*) through 
the diffeomorphism 

T,: Gxgxg* J-\r^) 

(x,y,0 ^ (x,F,Ad:r/ + ad;.e,0- 
Consider the G^-action on G x g x g* given, for each g E Gn and {x,Y, ^) E G x g x g* , 
hj g- {x,Y,^) = {gx, Y,^) and consider also the corresponding orbit space G x gx g* /Grj. 

Lemma 3 The diffeomorphism / f^) induces a new diffeomorphism 

T, : (G X g X g*) /G, ^ (G x g x g* x g*)^, (27) 

which maps the G^-orbit of the element {x,Y,^)EGxgx g* into the Gj^-orbit of the 
element (x, F, Ad*?^ + ady.^, .^) E J'^iji). 
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Proof. We have just to prove that is equivariant for the G^-action on J ^(77) and 
the G^-action on G x g x g* described above. Indeed, if (7 G G^ and {x,Y,^) E Gx qxq*, 



then T^{g ■ = {gx,Y, Ad*g^r] + ady^,0 = {9X,Y, Adl{Ad*gr]) + ady^,0 = 

{gx, Y, Adlv + ad^-e, = 09(T,(x, Y, 0)- □ 



Proposition 3 The point reduced space (G x g x g* x g*)^ is diffeomorphic to the space 
Orj X Q X Q* , where On denotes the coadjoint orbit of the element rj. 

Proof. It is clear that the map : (G x g x g*) /G^ — x g x g*, which takes a 
G,,-orbit of an element {x,Y,^) E Gxgx g* to a current point (Ad*?7, Y,^) E OrjXgxQ*, 



which gives us the diffeomorphism. □ 

The result from the previous proposition allows us to conclude now that (frj is a 
symplectomorphism, being 0^ = the symplectic structure on x g x g*. 

To finish, it is useful to notice that the map y?^ := o tt^ : J~^{rj) — )■ x g x g* 
is such that 



Furthermore, tp^ is surjective since (p^ is bijective and the projection vr^ is surjective. 

Besides the reduction of the phase space, the symplectic point reduction theorem 
has a dynamic counterpart, which will be addressed in the next subsection. 

4-2. Reduction of the dynamics 

We proceed with the analysis of the reduction of the dynamics of the Hamiltonian system 
of cubic polynomials (G x g x g* x g*, f2, H) described in subsection 13.31 We first present 
the natural reduction of dynamics that comes from the symplectic point reduction 
theorem. Then we shall perform this reduction as a dynamics on a Hamiltonian system 
on {O^ X g X g*,n,,) in the context of the previous subsection, that is, using the 
diffeomorphism fl28l) . 

Consider the Hamiltonian H given by f[T71) and the associated Hamiltonian vector 
field Xh defined by (fT8|) . Notice that H is invariant under the G-action defined by (12T|) . 
The symplectic point reduction theorem allows us to conclude the following: 

The flow ft of the Hamiltonian vector field induces a flow ft on the reduced 
space (G X g X g* X g*)^ defined by vr^ o o = ff^ o vr^. The vector field 
generated by the flow ft is Hamiltonian with associated reduced Hamiltonian 
function if^ defined uniquely by o vr^ = H o i^. Furthermore, the vector 
fields Xh and Xh„ are 7r„-related. 




(28) 



^,{x, Y, Adlv + ad^e, = (Ad:.r/, Y, 0- 



(29) 
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The triple {{G x g x g* x g*)^, fi^, if^) is called the reduced Hamiltonian system. We 
are interested now in characterizing the corresponding system on [O^ x g x q*,Q^). 
Namely, we shall determine the expression of the reduced Hamiltonian vector field when 
interpreted as a vector field on x g x g* regarding the description given at the end 
of the previous subsection. To effect this, one follows the steps below. 

From now on, wherever there is no confusion, we will denote an element Ad*?7 G 
hy 6, 6 := Ad*?]. Introduce the Hamiltonian function h := o (^~^ : x g x g* — )■ M, 
with defined by (1251) . 

Lemma 4 For each {6, Y, ^) E Or, x g x g* , we have 

h{e,Y,o = 0{Y) + lax^)- (30) 

Proof. Since the function ip defined by ^2^) is surjective, we know that the element 
(x, F, ^ + ady^, G J-^{r]) is such that iprj{x,Y,e + ad*Y^,^) = {9,Y,^). Thus, 
h{6,Y,C,) = (Ht) o '^ti){x,Y,6 + ady^,^). But from the definition of Hri, we have 
o TTri = H o i^^, SO usc [Tl\ ) and the result follows. □ 

Since is a symplectomorphism, we see that the Hamiltonian vector field 
associated to h is such that o = T<^^ o Xh^. Then, Xh o o vr^ = T<^^ o Xh^ o vr^, 
that is, Xh o (fri = T^prj o Xh^ o tt^, where ^pri is the function defined by (l29l) . Now, use 
the fact that Xh and X^^ are vr^-related, that is, T-k^, o X^ o = X^^ o tt^, to obtain 

XhO iprj = TiprjO Xh oir,. (31) 

We shall develop the expression of Xh for each point in (9^ x g x g*, using the relation 
now obtained and after we present a useful remark. 

Lemma 5 If{a,b,c,d) G T(^,y,e+adj,e,0'^"^('7); ihen 

T(x,Y,0+ad^C,OV5»?(«: b, c, d) = (&dT^L^_^J, b, d^ . (32) 

Proof. First we would like to clarify that the choice of the element (a, b, c, d) is related 
to the fact that T(^x,Y,e+iid*yi,0'J~^(^) ^ subset of T(^,Y,6»+adj,c,o x g x g* x g*) = 
T^G X g X g* X g*. Now, let (3 = (/3i, /32, Ad^^r/ + ad^/s,^) be a curve in J^\r]) 
satisfying the initial conditions /9(0) = {x,Y,9 + ady.^,,^) and /3(0) = {a,b,c,d). 
Then, we know that T'(^.,y,6i+ad^5,5)V^»;(o, 6, c, d) = d(v9^ o /3)(0)/dt, which is equal to 
^^^T^i(o)i^,(o)-i/3i(o)^^^i(o)^'/^2(0),/33(0)), that is, {adT^^^.^a^^b, d) as we wanted to 
prove. □ 

We are now able to prove the following result. 

Proposition 4 For each point [6, F, .^) G x g x g* , the dynamical vector field Xh of 
the Hamiltonian system {Orj x gx g*,r2^,/i) is given by 

Xh{e,Y,0 = {^^d*y9,X^,-9). (33) 
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Proof. We begin by noticing that due to the relations fl29p and fl3ip . we 
get Xh{9,Y,^) = {Tipn o XH){x,Y,d + ady^,0- Now according ([18]) we have 
Xh{x,Y,9 + ady^,^) = (TeLj,.!^, Xg, 0, — 6'). It remains only to show that 
T[x,Yfi+aAYi,i)VriiTeLxY, X^,Q, —6) = (ady 6^, X^, — 6') , but this follows from (132|) taking 
(a,6,c,d) = {TeLxY,X^,0,-e). □ 

Thus, the equations of Hamilton on the reduced manifold Orj x Q x g* are given by 



Remark 2 In remark [H we have proved the equivalence between the solutions of the 
equations of Hamilton flU) and the Euler- Lagrange equations f2^) . It is obvious that 
the reduced dynamics described by l[3^) is also related with the variational approach of 
cubic polynomials. In fact, an integral curve of the reduced Hamiltonian vector field (33) 
give rise to a curve that satisfies the second equation of the Euler-Lagrange system [2D^) . 
Indeed, writing the first equation of (^3^ as an equation on the Lie algebra (see the end 
of subsection \2. 1\ for details on notation), we get Xg + [Y,Xg] = 0. But, by the other two 
equations of [34\ ), we know that Xq = —Y. We conclude that a solution of the reduced 
system gives a solution of 



4-3. Invariants along the extremal trajectories 

Integrals of motion of a dynamical system are quantities that are conserved along the 
flow of that system and can be sometimes associated to symmetries of the system. 
A classical result due to Liouville says that a dynamical system on a phase space of 
dimension 2N is completely integrable if it admits N (almost everywhere) functionally 
independent first integrals in involution (i.e., their Poisson brackets all vanish). However, 
these situations are rather rare. In practice, one often deals with Hamiltonian systems 
which admits a non-abelian group of symmetries or an abelian group of symmetries in 
number less than the required to have complete integrability. If some special conditions 
are satisfied, the non-abelian set of independent integrals can lead us to the integrability 
of the system (theorem of Mishchenko-Fomenko). But, in most cases one naturally 
expects to find only a number of independent Poisson commuting invariants less than 
N, which will allow us to partially reduce the original system. 

The problem we are concerned with in this subsection is the preliminar analysis of 
the symmetries of the Hamiltonian system (O^ x fl x 0*, ^77, h), so that we can study the 
integrability of the system in a forthcoming work. We shall find, by using the classical 
Lie-Cartan theorem, a maximum number of functionally independent first integrals in 
involution. In the rest of this paper, for the sake of simplicity, we shall use the following 
notation: 




(34) 



Y+[Y,Y] = 0. 



(35) 



dimg = n 



and 



dima 



V ~ 



2m 



(36) 
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(recall that the dimension of the coadjoint orbit is always even), where obviously 2m < n. 
So, the dimension of the phase space of our Hamiltonian system shall be 

dim(C^ X X 0*) = 2(ra + m). (37) 

A function / : (9^ x g x — )■ M is an integral of motion of our dynamical system (with 
associated vector field X^) if [Xh{w)]{f) = 0, that is, [{df){w)]{Xh{w)) = 0, for all 
w e Or, X g X g* . It is important to notice that d/ : x g x g* — )■ T*{Or, x g x g*) 
is such that df{w) G T^Or, x g* x g C g x g* x g. In that sense, we shall assume the 
notation df{w) = {df{w)/d9,df{w)/dY,df{w)/di). 

The Hamiltonian function is naturally an integral of motion of the Hamiltonian 
system. So, the function ( l30l) . that is, 

k^h = e{Y) + \i [x^) 

is an integral of motion. But besides that, we are able to prove the following result: 
Proposition 5 The functions /j+i : (9^ x g x g* — )■ M defined by 

/j+i = {9 + ady^)(/li), with Ai a fixed basis element ofg, (38) 
are integrals of motion of the Hamiltonian system {On x g x g*, H^, h). 

Proof. Consider w = {9,Y,^) e O^ x g x g* C g* x g x g*, with 9 = Ad*?7, for some 
X E G. An elementary computation gives 

dli+i{w) = {Ai, -ad^^.^, adyAi) . (39) 

Then, by ([33]), we get [Xh{li+i)]{w) = ^(adyA^) - ^{adA^ - ^(adyA,) = 
{[X^, X^], Ai) = 0, which shows that [{dli+i){w)]{Xh{w)) = 0, for each i = 
proving that the given functions are invariant. □ 

Remark 3 Recall that in the context of the variational approach two invariants are 
known, ^ and ([^j- These invariants are related with the invariants now obtained. 
Indeed, it is simple to prove (using theoremU\) that li = Ii and 2/2 — X]r=i ^f+i = ^(^e)- 

It is immediate to see that any linear combination of the n + 1 integrals of motion 
described above is also an integral of motion of the system. A natural question arises: 
to extract, from the set of invariant functions, a maximal set of independent commuting 
invariant functions. In order to answer to this problem, we begin with some comments: 

- The independence between the invariants /j+i, i = l,...,n, defined by (!38|) is 
obvious, but these invariants are not in general in involution with each other. 

- On the other hand, it is clear that li commutes with each one of these functions 
/j+i, since /i coincides with the Hamiltonian function. 

Furthermore, we have the following two results: 

Lemma 6 Consider the invariants /j+i, for i = l,...,n, defined by f38\) . The 
Hamiltonian li is functionally independent of all the invariants li^i . 
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Proof. Consider the coordinate expression for the invariants, with respect to the 
natural basis taken from the orthonormal basis {Aj}i=i,...,n of the Lie algebra g: 

h = e, + Elk=i c}iy,i 



k 



Notice that for i > m, 6i is a function of 6i, 6m- We shall prove that d/i Ad/j+i 7^ 0, for 
a dense domain of (9^ x g x g* [i = 1, ...,n). But it is immediate to isolate the coefficients 
corresponding to dy^ A d9j, which are equal to 9jdy^ A d9j. Thus, the condition of the 
wedge product being equal to zero implies that 6j = 0, Vj, which is impossible for any 
regular value of 77 different from zero. 

If ?7 = it is immediate to check that the invariants are also independent, since in 
this situation the Hamiltonian does not depend on Y while the other invariants do. □ 

Lemma 7 Considering the Poisson structure on Oj^ x q x q* , the set {/i+i}i=i,...,n is 
endowed with a Lie algebra structure that makes it isomorphic to the Lie algebra g. 

Proof. Consider the orthonormal basis {Ai, An} of the Lie algebra g and represent 
by the structure constants of this Lie algebra for this basis. If w = {9,Y,^) G 
Or, X Q X g*, we know that {li+i,lj+i}{w) = [dli+i{w)] [Xi^^^{w)), with Xi^^-^ denoting 
the Hamiltonian vector field associated to Ij+i- 

In order to proceed with the proof, let us find now the expression of the Hamiltonian 
vector field ^z^+i of /j+i in a similar way to the one used in proposition H] for X^. To 
do this, first consider the function Lj^i :G'xgxg*xg*— t-R uniquely characterized 
by the identity Lj+i o = Z^+i o y?^, where y?^ is the surjective function defined by ( l29l) . 
More precisely, we have Lj+i{z) = ^! Ole=^_ad* 5 ~ /^(^i) where the element z is 

such that w = ^Pr,{z), that is, z = {x,Y,6 + ady^,.^) G J~^{ri). Then, the Hamiltonian 
vector field X^.^^ associated to the function Lj+i is related with Xi.^^ as follows 

X,^^,{w) = T^^,{Xl,,M) = fad^.L O^Xl^^^Xt^ , 

where we use lemma [5] choosing {a,b,c,d) = Xi^+^^iz) = (^1^+^, ^1^.+,, ^1^.+,, 
In order to completely determine we must to find the expression of the 

components Xl,^_^, ^ij+i "^^j+i' ^ computation analogous to the one done for 
the components of the vector field Xh in proposition [2] shows that Xl.^_^ = TgLxAj, 
-^Lj+i — adyy4j and ^tj+i — ^^*Aj^^ hence 

Xi^^^iw) = (ad:^/, adyA,, ad^i^e) • (40) 
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n 



i,j,k=l 

The structure constants of the Lie algebra generated by the functions i = 1, ...,n 
and the structure constants of q coincide, so the algebras are isomorphic. □ 

Let us summarize the situation: 

- We have n + 1 smooth independent functions, the integrals of motion /i, Z2, ^n+i, 
whose differentials are linearly independent on (9^ x g x g* (see lemma [6]). 

- The linear span of these functions is closed with respect to the Poisson bracket (see 
lemma [7] and recall that li is in involution with all the other functions). 

Thus, the linear span £ of these n + 1 functions has a structure of a finite-dimensional 
real Lie algebra, with dim£ = n + 1. This algebra is called the algebra of integrals. 
Finally, we shall use the S. Lie - E. Cartan (see [6]): 

Theorem 2 (S. Lie - E. Cartan) Consider a Hamiltonian system {M,u,H) with 
first integrals Fi,...,Fk such that {Fi,Fj} = aij{Fi, F^) . Let F : M ^ Ml' he the 
natural mapping generated by these set of integrals. 

Suppose that the point c & is not a critical value of the mapping F and that 
in its neighborhood the rank of matrix [aij) is constant. Then in a small neighborhood 
U C M'^' of c one can find k independent functions ipj : U ^ M. such that the functions 
(f)j = ipj o F : N ^ M, where N = F~^{U), satisfy the relations 



whereas the remaining brackets vanish. The number 2q equals the rank of the 

matrix [aij]. 

We are interested in the following consequence of the above theorem (see [6]): 

Remark 4 Under the hypotheses of Lie- Cartan theorem and using the notation above, 
there are k — q independent integrals in involution: 02, 04, 02g-2, 02g, 02g+i, 0fc- 
As a consequence, the original Hamiltonian system can be reduced, by the method of 
Poincare, to a system with minus k — q degrees of freedom than the original one. 

We consider now an open dense subset D in x g x g* where the functional 
independence of the n + 1 integrals is satisfied and where the skew-symmetric Poisson 
bracket matrix ({Zj,/^}), 1 < j < n + 1, has maximal rank. Notice that, by lemma [71 
the maximum rank of this Poisson bracket matrix coincides with the maximum rank of 
the matrix Mg(a) = (Mjj(a)), with Mjj(a) = Yllt=i^ij'^k^ for a = (ai,...,a„) G M" and 



{01,02} 



{02g-l,02j = 1, 
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where C^^ are structure constants of the Lie algebra g. We fix the notation (note that 
this rank is always even): 

rg := ^max<jgKnrankMg(a). (41) 

In consideration of remark |H the 2(m + n)-dimensional Hamiltonian system 
{Orf X X 0*, Clrf, h) admits n + 1 — rg functions defined on an open subset of D, which 
form a set of independent integrals of the motion in involution. Thus, we can expect 
the system to be reduced to a system of dimension equal to 2(m + rg — 1). 

Example 1 Consider the problem of cubic polynomials on the Lie group S0{3), which 
can be illustrated by the dynamic optimal control problem of the spherical free rigid 
body study by the authors of this paper in J^. In this case, it is well known that the 
coadjoint orbit 0,j corresponds to a 2-dimensional sphere with radius \\ri\\. (For the 
singular case t] = the orbit reduces to one point.) So, considering the non-singular 
case, the symplectic reduced manifold On x so (3) x so* (3) has dimension equal to 8 and 
rg = 1. Applying the above theorem, the corresponding reduced Hamiltonian system has 
3 independent invariants in involution and it can be reduced to a 2-dimensional system. 

As mentioned at the beginning of this section, these results on integrals of the 
motion may have important implications from the point of view of the integrability of 
the corresponding dynamical systems. Thus, they are relevant at the level of determining 
the cubic polynomials. In conclusion, by using the Lie-Cartan theorem, we are able to 
reduce the Hamiltonian system to a system with at least, n + 1 — rg degrees of freedom 
less. We note that if m + rg = 1, the system will be completely integrable. But, this 
condition happens only in trivial cases, when the algebra g is abelian or the coadjoint 
orbit is reduced to a point. 

We hope to address this integrability subject with detail in a future paper. 
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